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ABSTRACT 

We investigate the long-term variability exhibited by the X-ray point sources in the 
starburst galaxy M82. By combining 9 Chandra observations taken between 1999 and 
2007, we detect 58 X-ray point sources within the Z?25 isophote of M82 down to a 
luminosity of ^ 10'^'' ergs s""'^. Of these 58 sources, we identify 3 supernova remnant 
candidates and one supersoft source. Twenty-six sources in M82 exhibit long-term (i.e., 
days to years) flux variability and 3 show long-term spectral variability. Furthermore, 
we classify 26 sources as variables and 10 as persistent sources. Among the total 26 
variables, 17 varied by a flux ratio of > 3 and 6 are transient candidates. By comparing 
with other nearby galaxies, M82 shows extremely strong long-term X-ray variability 
that 47% of the X-ray sources are variables with a flux ratio of > 3. The strong X-ray 
variability of M82 suggests that the population is dominated by X-ray binaries. 
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1 INTRODUCTION 

Based on previous studies of X-ray sources in external galax- 
ies (see Fabbiano and White 2006; Fabbiano 2006), the X- 
ray populations with luminosity above typical Chandra de- 
tection limit (~ lO^'^ ergss"'^) of nearby galaxies are dom- 
inated by X-ray binaries (XRBs), which consist of an ac- 
creting compact object (neutron star or black hole) and a 
stellar companion. A few young supernova remnants (SNRs) 
may also be expected. X-ray binaries with low mass stellar 
companions are called low-mass X-ray binaries (LMXBs), 
which are usually associated with old stellar populations 
(bulges, globular clusters and early- type galaxies). On the 
other hand, high-mass X-ray binaries (HMXBs) with high 
mass companions are associated with young stellar pop- 
ulations like spiral arms and star-forming regions. More- 
over, some ofF-nuclear ultraluminous X-ray sources (ULXs) 
with Lx > 10''® ergs s~^ (Eddington luminosity of stellar 
mass objects) were found in several galaxies, especially star- 
forming galaxies. The nature of ULXs is still intensely de- 
bated. For star-forming galaxies, the source populations are 
also dominated by XRBs, especially HMXBs (Helfand and 
Moran 2001). Because of the association between HMXBs 
and short-lived young stellar populations, the amount of 
HMXBs may constitute a tracer of recent star formation. 

Today's X-ray source population studies in nearby 
galaxies are mostly based on X-ray luminosity functions 
(XLFs) and the intensity and spectral properties of indi- 
vidual sources. Furthermore, multiwavelength observations 



and studies of counterparts also provide important informa- 
tion. X-ray observations of a given galaxy at multiple epochs 
have provided opportunities to study long-term variabili- 
ties of X-ray sources in galaxies. Until now, various kinds 
of long-term intensity and spectral variability of different 
types of X-ray sources have been detected in both our own 
Galaxy and external galaxies. For example, black hole XRBs 
in our Galaxy are observed in five distinct spectral/temporal 
states, which are believed to be related to different accretion 
rates (Esin et al. 1997). Strong flux and spectral variability 
are shown during state transitions on long timescale. Com- 
paring variability properties of extragalactic X-ray sources 
with the well studied Galactic X-ray sources can assist in 
probing the nature of the sources, identifying the sources 
and understanding the source populations. 

In this paper, we report the long-term variability of 
the X-ray sources in M82 by using 9 archival Chandra data 
sets. M82, a prototypical star-forming galaxy with a dis- 
tance of 3.6 Mpc (Freedman et al. 1994), exhibits intense 
star-formation activity and hosts the most luminous ULX, 
M82 X-1 (Lx > 10*^ ergss-i). Although M82 is an impor- 
tant nearby starburst galaxy, not much work has been done 
on the X-ray source population and the long-term variability 
of the whole galaxy due to its crowded source distribution, 
a strong absorption (A'h ~ 3 x 10^^ cm"'^) and the intense 
and non-uniform diffuse emission in the central region of 
M82. X-ray studies of M82 mostly focused on M82 X-1, the 
prime candidate of intermediate-mass black hole (IMBH). 
For the X-ray source population, Zezas et al. (2001) de- 



© 2010 RAS 



2 Yi-Kuan Chiang and Albert K. H. Kong 



Table 1. Chandra Observation Log 



No. 


Obs. ID 


Date 


MJD-51300 


Exp. 


Instrument 


1 


361 


1999/9/20 


141.47 


33.7ks 


ACIS-I 


2 


1302 


1999/9/20 


141.88 


15.7ks 


ACIS-I 


3 


1411-1 


1999/10/28 


179.18 


36.3ks 


HRC-I 


4 


1411-2 


2000/1/20 


263.64 


17.8ks 


HRC-I 


5 


2933 


2002/6/18 


1143.78 


18.3ks 


ACIS-S 


6 


5644 


2005/8/17 


2299.04 


75. Iks 


ACIS-S 


7 


6361 


2005/8/18 


2300.66 


19.2ks 


ACIS-S 


8 


8189 


2007/1/9 


2809.34 


61.6ks 


HRC-S 


9 


8505 


2007/1/12 


2812.09 


83.6ks 


HRC-S 



NOTE. — Obs. ID 1411 has two observations merged in one 
event list. We used a time filter to separate the observations. 

tected 24 point sources in M82 and 12 of them varied by 
factors of 0.2 to 7 on timescales of 1 to 6 months. An XLF 
of the sources was then constructed and can be fitted with 
a power law with a slope of ^ —0.4. Kilgard et al. (2002) 
studied the X-ray point sources of three starburst galax- 
ies and four nonstarburst spiral galaxies. They presented 
an XLF of more than 30 sources detected in M82 with a 
slope of ~ —0.5 and further showed that the XLFs of star- 
burst galaxies tend to be flatter than those of spiral galax- 
ies. Griffiths et al. (2000) studied the hot plasma and the 
point sources in M82. They resolved about 20 compact X- 
ray sources and suggested that many of which could be black 
hole HMXBs. Now, much more X-ray observations of M82 
at multiple epochs have been done, allowing us to study the 
X-ray population and the long-term variability in detail. 

In § 2 we describe the Chandra observations. We present 
the data analysis in § 3 and the results in § 4. A discussion 
of the results is given in § 5. 



2 OBSERVATIONS 

Since the X-ray sources in the central region of M82 are 
crowded, a high spatial resolution instrument like the Chan- 
dra X-ray Observatory is needed to resolve the images. M82 
has been observed 13 times between the year of 1999 and 
2007 with Chandra; some of them were designed to observe 
M82 X-1 with special configurations. To study the long-term 
variability of the whole galaxy, we used 9 Chandra data sets 
with relatively long exposure time and large field-of-view 
from the Chandra Data Archive. The details of these obser- 
vations are given in Table 1. Among these 9 observations, 
observation 3, 4, 8 and 9 were performed with the High 
Resolution Camera (HRC-I or HRC-S). The rest were per- 
formed with the Advanced CCD Imaging Spectrometer ar- 
ray (ACIS-I or ACIS-S). For Obs. ID 1411, there are two 
separate observations within the same event list. Therefore, 
we used a time-filter to split them (observation 3 and 4) and 
analyzed the data separately. 



3 DATA ANALYSIS 

The data reduction and analysis were done with CIAO, Ver- 
sion 4.2 (CALDB Version 4.3.0) and HEAsoft, Version 6.3. 
The ACIS and HRC data sets were processed with the stan- 
dard procedure by the Chandra X-ray Center (CXC). We 
obtained the fully processed science products (level 2 event 



files) of these observations. For the ACIS data sets, we only 
used events between 0.3 — 7 keV. In order to remove the 
differences in the coordinates of these Chandra data sets, we 
registered them using 3 bright X-ray point sources. 

3.1 Source detection 

With WAVDETECT, a wavelet detection algorithm imple- 
mented within CIAO, we can search for and determine the 
position of the point sources. To process the images of the 
crowded central region, the 0.49 arcsec pixel size of ACIS 
images is insufficient. In contrast, HRC images have smaller 
pixels (0.13 arcsec) which oversample the PSF of an on axis 
point source. Therefore, HRC has the highest spatial reso- 
lution of Chandra. For all the ACIS data sets, we generated 
the subpixel images with a pixel size of 0.123 arcsec (=1/4 
ACIS pixel) to reach the resolution of HRC images. 

We applied WAVDETECT on aU the 0.3 - 7 keV 1/4 
ACIS subpixel images and the HRC images at five scales (1, 
2, 4, 8 and 16 pixels). We set the detection thresholds to be 
the inverse of the total pixel number of the images, which 
is equivalent to stating that the expected number of false 
sources per field is one. However, there are some spurious 
detections due to the strong diffuse emission in the central 
region of M82. After visual inspection, we obtained a can- 
didate source list by combining the WAVDETECT source 
lists generated for each observation. We used circular source 
extraction regions of varying radii between 1" to 4" to cover 
more than 90% energy of the sources. For each source, a 
nearby source free region was extracted as background. The 
background extractions need to be treated especially care- 
fully in the central region since the diffuse emission is not 
uniform. After intensity and spectral analyses, source can- 
didates with a flux below Icr were first rejected. A source 
candidate embedded in a relatively irregular part of the dif- 
fuse emission will be further considered to be a background 
contamination if its color is as soft as the diffuse hot gas. In 
general, sources within the optical _D25 isophote (~ 11.6' x 
5.6'; ~ 12 kpc x 6 kpc) are thought to be associated with 
M82. The background contribution was further estimated. 

3.2 Light curves 

To estimate the X-ray flux of each source, we assume an 
absorbed power-law model with a photon index of 1.7 and 
a column density of hydrogen A^h = 3 x 10^^ cm~^ (the 
typical value we obtained from the spectral flttings of the 
bright sources in M82). We generated the redistribution ma- 
trix files (RMFs) and the auxiliary response files (ARFs) of 
the sources in each observation. Applying the response files 
and the spectral model described above, we obtained the en- 
ergy fluxes in 0.3 — 7 keV using a tool called modelflux. We 
further converted these absorbed fluxes into intrinsic fiuxes 
with PIMMS (Portable, Interactive Multi-Mission Simula- 
tor). A distance of 3.6 Mpc to M82 (Freedman et al. 1994) is 
adopted for the estimation of luminosities. If the net counts 
of a source in an observation is below Icr, we obtained the 
90% upper limit of the counts (Gehrels 1986) and the corre- 
sponding fiux. To deal with the Poisson and binomial statis- 
tics inherent to small numbers of counts, all the errors of 
counts are calculated using Gehrels' error approximation 
(Gehrels 1986). 
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3.3 Hardness ratios 

Since many of the sources are too dim, it is inappropriate 
to fit the spectra of these faint sources with small degrees 
of freedom. To obtain a rough picture of spectral properties, 
we computed the hardness ratios of all the sources in each 
ACTS data set (observation 1, 2, 5, 6 and 7). The soft color 
is defined as 

(1) 



HRl = 



T 



and the hard color as 
H-M 



HR2 



T 



(2) 



where S, M and H are the net counts in the soft (0.3 — 1 
keV), medium (1 — 2 keV), and hard (2 — 7 keV) bands, 
and T are the net counts in the whole energy band(0.3 — 7 
keV) . The errors of hardness ratios were derived from error 
propagation. 



3.4 Variability 

To estimate long timescale flux variability, we calculated the 
maximum significance of the difference in fluxes between any 
two observations. 



max 



(3) 



where Fi is the energy flux of a source in the ith observa- 
tion and is the corresponding error. The 90% upper limit 
was derived when a source is below la detection. Snux rep- 
resents the confidence level of fiux variation. In this paper, 
variations with 3a significance are thought to be detected. 
However, it is inevitable that the significance parameter of 
a source is affected by the fiux level it reaches, since the 
fractional error of flux of a faint source is larger than that of 
a brighter source. In other words, variations of faint sources 
are hard to be detected even if these sources went through 
intrinsic variations between observations. 

To quantify the intensity of variations, we also com- 
puted the ratio of fluxes. 



R 



F ■ ' 

-L mm 



(4) 



where -Fmax and -Fmin are the maximum and minimum fluxes 
of a source among these 9 observations. R could be a number 
with a standard error if the lowest flux is detected or a 90% 
lower bound if the lowest flux is an upper limit. 

In addition, we checked intra-observation variability for 
the sources in all the data sets via the Gregory-Loredo al- 
gorithm (Gregory and Loredo 1992). Source variability is 
evaluated from the arrival times of photon events with m 
bins per observation, where m runs from 2 to mmax. Frac- 
tional areas considering the effects of dither motion and bad 
pixels were corrected. Based on the Gregory-Loredo prob- 
ability, the odds ratio and the deviation of the light curve 
from the mean value, an integer variability index could be 
assigned for the source, which provides a measurement of 
intra-observation variability in various time scales. 

For long-term spectral variability, we computed the sig- 
nificance of the difference in soft color as 



|H R1, - HRIj I 



•Shri = niax 



and the significance of the difference in hard color as 
|HR2, -HR2, i 

OHR2 — max 



(5) 



(6) 



It should be noted that we can obtain spectral information 
only in 5 observations (observation 1, 2, 5, 6 and 7) which 
were performed with ACIS-I or ACIS-S. These 5 ACIS ob- 
servations span over 6 years. 



4 RESULTS 

4.1 Source detection 

After the processes we described in § 3, we detected 58 X-ray 
point sources within the optical 1)25 isophote of M82 in these 
9 Chandra observations during a span of 7 years. Depending 
on sources' location and exposure time, the detection limit 
is 8 X 10'''' ergs s~^ in the outer region and increases to 
3 X 10^^ ergss"^ in the central region due to the strong dif- 
fuse emission. Among the total 58 X-ray point sources we 
detected, 25 are in the central 1' x 1' region and 33 sources 
are in the outer region with fluxes down to a fainter limit. 
One of the 25 detections in the central 1' x 1' region and 
6 of the 33 detections in the outer region are expected to 
be background sources. The background contribution we es- 
timated is based on the results of the Chandra deep fleld 
north survey (Brandt et al. 2001), which provides an esti- 
mate of background sources (AGNs or starburst galaxies) 
density above a given flux limit. 

We also compared the positions of our sources with a 
list of 37 supernova remnants provided by the deep MERLIN 
5GHz radio observation (Fenech et al. 2008), which covers 
the central 1 kpc region of M82. Two of the X-ray sources 
(CXOU J095550.7-H694044 and CXOU J095552.8+694045) 
have corresponding radio-selected SNRs located within their 
positional errors. For the historical type II supernova SN 
2004am (Singer et al. 2004; Mattila et al. 2004; Beswick et 
al. 2004), there is an X-ray source CXOU J095546.7-f694038 
associated with the optical position. 



4.2 Variability 

As we described in § 3.2, the 0.3 — 7 keV unabsorbed fluxes 
were obtained by assuming an absorbed power-law model 
with a photon index of 1.7 and a column density A^h = 
3 X 10^^ cm~^. We further evaluated the variation of the 
derived unabsorbed fluxes with different spectral parame- 
ters. Varying photon index from 1.7 to 1.2 and 2.0 for the 
sources in observation 1 results in flux changes within 20%. 
Combined with the fact that most of the sources in M82 do 
not undergo spectral variation significantly (details will be 
presented in this section), the long-term variability of the 
sources is insensitive to the assumed spectral parameters. 

Based on the significance parameter of long-term flux 
variability 5flux deflned in § 3.4, we classify the X-ray point 
sources into persistent sources if Saux < 3 or variables if 
•Sflux > 3. For some faint sources with fluxes close to the 
detection limit, we can only obtain the lower limit of their 
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Table 2. Chandra source properties for M82 



No 




Source Name 
(CXOU J) 


R. A. Dccl. 
(J2000) 


Max. Ly^^ 
(10^** orgss^^) 




R" 


Shri'' 


Sun 2" 


Notcs^ 


1 




095510.7-1-693955 


148.79455 


69.665369 


2.68 


2.78 


3.64±1.41 


0.58 


1.2 


P 


2 




095527.2-1-693923 


148.86338 


69.656503 


11.71 


5.59 


2.23±0.3 


1.5 


1.67 


V 


3 




095527.2-1-694050 


L48. 86345 


69.680558 


0.58 


>2.58 


>8.48 








4 




095529.1-1-694027 


L48. 87116 


69.674172 


0.88 


>2.68 


>6.45 


1.24 


1.41 




5 




095529.2-1-694212 


148.87180 


69.703353 


1.14 


>1.15 


>3.09 


0.78 


0.35 




6 




095530.2-1-694239 


148.87570 


69.710836 


0.81 


>0.85 


>2.37 


0.07 


0.41 




7 




095530.9-1-694234 


148.87886 


69.709403 


1.15 


>1.68 


>5.22 


0.41 


0.4 




8 




095531.9-1-693958 


L48. 88280 


69.665981 


0.87 


>1.45 


>3.98 


0.43 


0.52 




9 




095532.7-1-694001 


148.88633 


69.667075 


2.98 


4.35 


21.5±14.93 


0.82 


0.29 


V3 


10 




095534.2-1-693943 


148.89275 


69.662022 


0.85 


>1.59 


>4.90 


0.33 


0.53 




11 




095534.5-1-693824 


148.89392 


69.639917 


0.95 


>1.51 


>2.77 


0.93 


0.33 




12 




095537.9-1-694057 


148.90792 


69.682597 


0.90 


>1.64 


>7.85 


0.56 


1.16 




13 




095538.0-1-694030 


148.90853 


69.674875 


6.02 


4.58 


2.16±0.32 


1.19 


1.64 


V 


14 




095541.0-1-693928 


L48. 92068 


69.657761 


2.76 


>4.43 


>25.41 






V3 


15 




095542.8-1-694033 


148.92850 


69.675839 


0.81 


>1.18 


>2.81 


0.4 


0.82 




16 




095544.1-1-693958 


148.93383 


69.666206 


0.61 


>1.18 


>2.74 


0.47 


0.23 




17 




095546.2-1-694026 


L48. 94247 


69.673997 


5.47 


2.79 


1.75±0.31 


2.32 


1.14 


P 


18 




095546.6-1-694041 


148.94403 


69.677994 


98.16 


37.57 


4.69±0.36 


3.56 


5.07 


V3, SV, IOV{6,7) 


19 




095546.7-1-694038 


148.94453 


69.677175 


5.72 


5.89 


2.34±1.03 


0.98 


0.86 


SN 2004am, V 


20 




095547.4-1-694036 


148.94768 


69.676733 


46.06 


>21.18 


>884.19 






T, IOV(3) 


21 




095547.5-1-694100 


148.94777 


69.683258 


4.31 


2.68 


1.45±0.27 


1.07 


1.72 


P 


22 




095548.8-1-694044 


148.95328 


69.678783 


0.84 


>0.61 


>1.64 


0.79 


0.02 


IOV(7) 


23 




095549.4-1-694044 


148.95594 


69.678783 


2.88 


1.41 


1.99±1.23 


1.68 


1.12 


P, IOV{8) 


24 




095550.0-1-693934 


148.95822 


69.659514 


4.34 


6.41 


5.31±1.39 


0.97 


1.66 


V3 


25 




095550.1-1-694046 


148.95889 


69.679569 


1581.69 


181.28 


28.25±0.82 


8.75" 3.47'' i 


.37' 4.76^ 


M82 X-1, V3, SV, IOV(6) 


26 




095550.3-1-694022 


148.95978 


69.672908 


8.36 


5.77 


2.84±0.51 


0.52 


1.58 


V 


27 




095550.4-1-694036 


148.95987 


69.676767 


11.34 


>15.84 


>78.22 


0.31 


0.87 


T 


28 




095550. 7-H694044 


148.96105 


69.678783 


20.62 


5.39 


1.75±0.25 


2.18 


3.12 


RSN, V, SV 


29 




095551.0-1-694045 


L48. 96243 


69.679194 


108.98 


98.6 


65.7±25.61 


1.15 


2.21 


T 


30 




095551.3-1-694044 


148.96361 


69.678853 


35.63 


12.28 


2.85±0.31 


3.17 


1.64 


V 


31 




095551.5-1-694036 


L48. 96450 


69.676667 


22.33 


5.79 


1.93±0.31 


0.93 


1.72 


V 


32 




095551.9-1-694050 


148.96607 


69.680458 


1.36 


>1.97 


>6.15 


1.91 


1.41 




33 




095552. 0-H694043 


148.96647 


69.678614 


4.90 


>4.93 


>3.99 


0.8 


1.17 


V3 


34 




095552.3-1-694054 


148.96774 


69.681586 


24.27 


>24.59 


>87.42 


1.12 


2.68 


T 


35 




095552.6-1-694047 


148.96932 


69.679808 


3.80 


2.85 


5.62±4.75 


1.36 


2.81 


P 


36 




095552.8-1-694045 


148.96981 


69.679297 


5.43 


2.81 


4.25±4.11 


1.73 


1.27 


RSN, P 


37 




095553.2-1-694049 


148.97168 


69.680219 


5.50 


4.11 


4.19±1.89 


1.58 


1.92 


V3 


38 




095553.4-1-694102 


148.97256 


69.683875 


2.20 


1.81 


2.18±1.16 


0.97 


1.02 


P 


39 




095553.4-1-693956 


148.97257 


69.665461 


2.54 


2.12 


3.35±3.16 


1.25 


0.8 


P, IOV{2) 


40 




095553.8-1-694050 


148.97414 


69.680664 


9.99 


>8 


>7.77 


1.73 


1.44 


V3 


41 




095554.2-1-694040 


148.97571 


69.677692 


10.45 


11.24 


11.91±5.26 


2.24 


3.81* 


V3 


42 




095554.7-1-694101 


148.97787 


69.683569 


10.65 


6.29 


2.19±0.27 


1.59 


1.41 


V 


43 




095555.1-1-694028 


148.97945 


69.674481 


1.23 


>1.14 


>3.07 


0.94 


0.24 




44 




095600.3-1-694046 


149.00139 


69.679400 


1.26 


>2.8 


>12.46 


0.32 


0.07 


IOV(9) 


45 




095601.3-1-694111 


149.00553 


69.686469 


18.30 


>17.74 


>112.17 






T 


46 




095601.7-1-694309 


149.00705 


69.719167 


1.11 


>1.63 


>2.94 


0.5 


0.75 




47 




095606.0-1-694031 


149.02500 


69.675328 


0.97 


>1.91 


>11.59 


0.66 


0.16 




48 




095608.1-1-694140 


L49. 03360 


69.694322 


5.94 


>8.84 


>126.36 






T 


49 




095614.8-1-694249 


149.06164 


69.713617 


1.10 


>3.72 


>9.16 


0.65 


0.48 


V3 


50 




095615.1-1-693952 


L49. 06283 


69.664347 


1.25 


>1.67 


>3.42 


0.26 


0.67 




51 




095615.2-1-694142 


L49. 06334 


69.695097 


1.09 


>4.01 


>13.71 






SSS, V3 


52 




095619.7-1-694146 


149.08225 


69.696147 


6.89 


3.71 


1.91±0.29 


0.76 


1.06 


V 


53 




095621.2-1-694225 


149.08821 


69.706939 


1.35 


1.69 


3.37±1.74 


0.43 


0.29 


P 


54 




095621.4-1-693901 


149.08920 


69.650153 


0.65 


>0.32 


>1.54 


0.21 


0.29 




55 




095622.7-1-694135 


149.09463 


69.692961 


1.13 


>1.85 


>14.26 


0.41 


0.01 




56 




095631.1-1-694220 


149.12977 


69.705611 


2.88 


>2.47 


>2.87 


1.01 


0.91 




57 




095633.8-1-694352 


149.14080 


69.731192 


1.60 


>1.65 


>2.92 


0.71 


0.56 




58 




095637.6-1-694137 


149.15686 


69.693731 


1.25 


1.34 


1.62±1.05 


0.87 


1.19 


P 


" The highest luminosity in 0.3 - 7 keV of a source among the 9 observations we used, assuming an absorbed power-law model 
with a photon index of 1.7 , A^h = 3 x 10^^ cm"^ and d=3.6 Mpc. 

Sfl^, = max,,j(|F, - Fj\/^a%.^ + af,^ ) 

R = i^max/i^min with a Standard error if the lowest flux is detected or a 90% lower bound if the lowest flux is an upper limit. 
Shri = inaxi^j ( HRli — HRlj |/ ^/'^'hri "^hri ) °^ ^ Chandra ACIS observations. 

S'hr2 = maxi^j (|HR2i ~ HR2j / ^cr^jj^j ''^hr2 ) of the 5 Chandra ACl^ observations. 

^ P: persistent source, V: variable, V3: variable with flux ratio > 3, T: transient candidate, SV: spectral variable, SSS: supersoft 
source candidate, RSN: source coincides with a radio-selected supernova remnant in the list of Fenech et al.(2008), IOV(Obs. 
No.): intra-observation variability detected in the specified observation. 

' The ^HRi & S'hr2 with a different extraction of colors in observation 5 using a pile-up spectral fitting. This source suffers 
from mild pile-up in observation 5 with a pile-up fraction of ~ 26%. 

^ The Shri ^ Shr2 excluding observation 5. 



^ Considering the strong contamination from the diffuse emission, this significant level is overestimated. 
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Sflux- Therefore, it would be impossible for us to classify a 
source if its lower limit of Saux < 3. To indicate the amount 
of variability, we also mark the variables whose flux ratio 
7? > 3. Furthermore, if a variable showed a strong variability 
with R > 50, we consider it to be a transient candidate. To 
describe the long-term variability of the whole galaxy, we use 
the fraction of variables with i? > 3 (excluding unclassified 
sources) to be an indicator. 

Among the total 58 x-ray point sources in M82, there 
are 10 persistent sources, 26 variables and 22 unclassified 
sources. Seventeen of the 26 variables showed a variability 
with _R > 3. Six of the 26 variables are transient candidates. 
Excluding the unclassified sources, 47% of the detected X- 
ray sources in M82 are variables with _R > 3 on timescale 
of months-years. The long-term variability properties of the 
sources in M82 are given in Table 2. 

For the central 1' x 1' region of M82, the population 
of long-term flux variability is similar to that of the whole 
galaxy. 

Intra-observation variability plays a limited role in the 
detected fiux variances between observations. Sources with 
Gregory-Loredo variability index less than 3 are considered 
to be persistent within observation. Of all the sources in 
M82, only 6 and 1 showed short-term variability in 1 and 2 
observations respectively. We note these sources with corre- 
sponding observations in Table 2. 

Furthermore, we estimated the spectral variability us- 
ing the significance of the difference in hardness ratios. We 
consider a source to be a spectral variable if Shri > 3 or 
5'hr2 > 3. Originally, we found 5 sources satisfied these crite- 
ria but 2 spurious cases have been dismissed. Spurious vari- 
ation could be detected in a transient like source embedded 
in a region of strong and irregular diffuse emission. Although 
the high error of flux when a source was faint should sup- 
press the confldence of its spectral variation, contamination 
from diffuse emission could be relatively serious in its faint 
states. Then, a bright outburst could make the composite 
error small compared to the color difference and a spurious 
detection may be selected. If a source suffered from strong 
contamination of the diffuse emission, we expect that the 
error should be greater than the standard error, and the de- 
tection slightly over our criteria should be dismissed. After 
all, 3 sources exhibited long-term spectral variability with 
SwRi > 3 or 5hr2 > 3. All of these 3 spectral variables also 
show variability in their X-ray flux. 

Due to the limiting statistics, it is inevitable that the 
transient candidates and the spectral variables we found 
tend to be luminous. Sources with relatively low flux and 
sources within the region of strong diffuse emission hardly 
reach our criteria of transient candidates and spectral vari- 
ables due to this selection effect. 

4.3 Individual sources 

For the sources which are classifled as variables, transient 
candidates or spectral variables, various kinds of variability 
properties are shown. Some of them may be related to spec- 
tral/temporal state transitions of black hole XRBs. How- 
ever, it is usually hard to confirm because of the lack of low 
state spectra with our detection limit. X-ray power spec- 
tra and studies of radio counterparts are also very useful to 
probe the nature of the sources. 



Some of the transient candidates exhibited outburst 
once during these 9 Chandra observations. One of the tran- 
sients showed a lightcurve which is consistent with a fast rise 
followed by an exponential decay (FRED). Two transients 
showed recurrent outbursts on long timescale. In addition to 
the 6 transient candidates, some of our sources also showed 
transient-like behavior but did not meet our criteria for be- 
ing classified as a transient candidate. We also found a super- 
soft source (SSS), which only appeared in one observation 
with all its counts in the soft band and was below detection 
limit in the other observations. For the 3 spectral variables, 
we constructed the color-luminosity diagrams. The detailed 
properties of the individual sources will be presented in this 
section. 

4-3.1 Transients 

CXOU J095547.4+694036, an X-ray transient candidate 
which showed a long-term lightcurve consistent with a fast 
rise followed by an exponential decay (FRED), a typical 
profile of X-ray novae. This source reached its peak lumi- 
nosity (4.6 X 10"^^ ergss"^) in observation 3 and followed by 
a lower luminosity of 1.3 x 10^^ ergss"^ three months later. 
The detailed profile of its flux decay is unclear due to in- 
sufficient data points. We do not have spectral information 
because the 2 observations which caught the outburst were 
performed with HRC-I. 

CXOU J095552.3-f694054, an X-ray transient candi- 
date which was luminous in the first two observations with 
Lx > 10^^ ergs s~^. And then it gradually decayed to the 
detection limit on time scale of years. It is not clear whether 
it has a fast rise profile or not. The spectra of CXOU 
J095552. 3+694054 in observation 1 can be well fitted with 
an absorbed power-law model with a photon index of 1.8. 

CXOU J095601. 3+694111, another transient candidate, 
originally showed luminosities below the detection limit. 
However, it appeared in the last two HRC observations with 
Lx > 10^® ergss"^ and no spectral information is available. 
The profile of the outburst is not clear. This source is not 
covered by the field-of-views of observation 6 and 7. 

CXOU J095608. 1+694140, an X-ray transient can- 
didate suddenly appeared in observation 5 with Lx > 
10''* ergs s~^ and seems to be stable in fiux hereafter. The 
X-ray colors of CXOU J095608. 1+694140 also did not vary 
after observation 5. 

The other two transient candidates we found in M82, 
CXOU J095550.4+694036 and CXOU J095551. 0+694045 
showed recurrent outbursts between 1999 and 2007. One of 
the most luminous X-ray sources near the center of M82, 
CXOU J095551. 0+694045, has been studied using roughly 
the same observations in the work of Kong et al. (2007). 
The X-ray luminosity of CXOU J095551.0+694045 reaches 
10*'' ergs s~^ in outbursts. When it turned into quiescence 
in observation 1, 2, 4, 8 and 9, basically it cannot be 
detected by WAVDETECT. The detected luminosities (~ 
10'^* ergs s~^) in observation 1, 2, 4, 8 and 9 that we 
present in the lightcurve (Figure 1) may due to contami- 
nation of the diffuse emission and nearby luminous point 
sources. Indeed, the detection limits for this position are 
just about 10^* ergs s~^. Our lightcurve of this ultralumi- 
nous X-ray transient is consistent with which shown in Kong 
et al. (2007). For the other two recurrent X-ray transients. 
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Figure 1. Long-term X-ray lightcurve of the transient candidates. The luminosities are determined by assuming an absorbed power-law 
model with a photon index of 1.7 , A^h = 3 x 10^^ cm~^ and d=3.6 Mpc. For non-detections, downward triangles represent their 90% 
upper limits. 



one of them (CXOU J095550.4-I-694036) has spectra avail- 
able in ACIS observations. The spectra can be fitted with 
an absorbed power-law with a photon index of ~ 3.7 and a 
reduced of 0.9, indicating that the source is in high/soft 
state. 

In addition to these 6 transient candidates, some 
dim sources may reach our criteria of transient can- 
didates if there was no diffuse emission. One case is 
CXOU J095553. 8+694050, a recurrent source located in the 
strong diffuse emission region. Another transient-like source. 



CXOU J095615. 2+694142, is the only SSS we found in M82. 
It only appeared in observation 5 with its total counts of 30 
all in the soft band. The colors and the transient nature of 
this source are very similar to that of typical SSSs in nearby 
galaxies (Di Stefano and Kong 2003). We also estimated 
the luminosity of the SSS by assuming a blackbody model 
with temperatures of 50 eV and 100 eV. The estimated lu- 
minosity strongly depends on the assumed column density 
of hydrogen. The accurate column density should be located 
between the value of the Milky Way absorption in this direc- 
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Figure 2. Hard color as a function of luminosity for sources that show significant variation in their hardness ratios and are classified as 
spectral variables. The numbers next to the data points indicate the Obs. No. 



tion (A^H = 4 X 10^° cm~^, Dickey and Lockman, 1990) and 
the typical absorption (A'h = 3 x 1Q^^ cm~^) obtained from 
the spectral fitting of the bright sources in the central region 
of M82. With the Galactic absorption (iVn = 4 x 10^" cm~^) 
and a blackbody model with a temperature of 50 eV, the ab- 
sorbed and unabsorbed 0.3-7 keV luminosities are 1.90 x 10'^^ 
and 3.90 X 10^^ ergs s~^, respectively. For a blackbody model 
with a temperature of 100 eV, the absorbed and unab- 
sorbed 0.3-7 keV luminosities are 1.27 x lO^'^ and 1.99 x lO^'^ 
ergss"^, respectively. If we assume a higher typical absorp- 
tion of the central bright sources in M82 (A^h = 3 x 10^^ 
cm~^), the expected unabsorbed 0.3-7 keV luminosity of 
this SSS could reach 7.68 x 10*^ ergss"^ using a black- 
body model with a temperature of 50 eV and 5.21 x 10*° 
ergss"^ for a black body temperature of 100 eV. If the ab- 
sorption is indeed much higher than the Galactic one, CXOU 
J095615.2-f 694142 should be classified as an ultraluminous 
SSS like the one in MlOl (Kong et al. 2004; Kong and Di 
Stefano 2005). 



servation 5. However, medium pile-up with a pile-up frac- 
tion of ~ 26% from which it suffered could still impair 
the quality of spectral analysis. Excluding observation 5, 
the significance parameter of the difference in soft color 
(Shri) reduces from 8.75 to 3.47. The significance param- 
eter of hard color (S'hr2) also reduces from 8.37 to 4.76. 
M82 X-1 still seems to undergo spectral variation during 
the other ACTS observations (Obs. 1, 2, 6 and 7). The 
color-luminosity diagram and the spectra reveal a transi- 
tion from a low/hard state to a high/soft state. Detailed 
analysis and discussion of the spectral properties of M82 
X-1 were presented by Agrawal and Misra (2006). Another 
spectral variable, CXOU J095550. 7-1-694044, showed a in- 
verse color-luminosity correlation to that of M82 X-1. In 
fact, CXOU J095550.7-f 694044 is a luminous young X-ray 
supernova remnant candidate near the center of M82, we will 
discuss it more in the next subsection. The other spectral 
variable we found in M82, CXOU J095546.6+694041, did 
not show a correlation between luminosity and color clearly. 



4.-3.2 Spectral variables 

M82 X-1, CXOU J095550. 1-^694046, showed strong X-ray 
spectral variability in the data set of observation 5 due to 
pile- up effect. We used the pile- up model to fit the spec- 
tra and extracted the hardness ratios for M82 X-1 in ob- 



4-3.3 Supernova remnants 

CXOU J095546. 7-1-694038, a source which is at the position 
of a historical supernova SN 2004am (Singer et al. 2004; 
Mattila et al. 2004; Beswick et al. 2004). SN 2004am was 
discovered in M82 on March 5, 2004 by Lick Observatory 
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Figure 3. Long-term lightcurve of CXOU J095546. 7+694038, an 
X-ray source at the position of a historical supernova SN 2004am. 
In the observations after the supernova explosion (MJD— 51300 ~ 
1650), the X-ray luminosity of it increased by a factor of ~ 50%. 
The luminosities are determined by assuming an absorbed power- 
law model with a photon index of 1.7 , Nn = 3 X lO^'^ cm"'^ and 
d=3.6 Mpc. 



Supernova Search (LOSS) and the explosion is believed to 
take place at the end of 2003. SN 2004am was determined 
as a type II supernova. Although the positions of the X- 
ray source CXOU J095546. 74-694038 and SN 2004am match 
well, this X-ray source has been detected before the super- 
nova explosion. It suggests that CXOU J095546. 7-1-694038 
is not the X-ray counterpart of SN 2004am. However, the 
lightcurve of CXOU J095546. 74-694038 (Figure 3) shows 
that in the observations after the supernova explosion, its 
X-ray luminosity increased by a factor of ~ 50%. It is not 
clear whether the increase of X-ray luminosity is due to the 
appearance of the supernova remnant of SN 2004am along 
the line of sight to the original X-ray source or not. We also 
check the Chandra spectra before and after the supernova 
event, the quality of these spectra is not good enough for us 
to explore the nature of the source. 

CXOU J095550.7-F694044, an ultraluminous SNR can- 
didate near the center of M82, has a strong radio coun- 
terpart known as 41.95-1-575 (McDonald et al. 2001, 2002; 
Fenech et al. 2008). It has been suggested to be a young 
supernova event taking place within a high density molecu- 
lar cloud (McDonald et al. 2001). The X-ray spectral prop- 
erty of CXOU J095550. 7+694044 has also been discussed 
in the work of Kong et al. (2007). In the archival Chan- 
dra ACIS data sets we used (observation 1, 2, 5, 6 and 7), 
it showed a spectral variation in HR2. The X-ray color of 
CXOU J095550. 7+694044 seems softer in the observations 
when CXOU J095550. 7+694044 is relatively dim. 

CXOU J095552.8+694045 is also located within the er- 
ror circle of a radio-selected supernova remnant in the cata- 
log of the deep MERLIN 5GHz radio observation (Fenech et 
al. 2008). The flux of this X-ray source is close to the detec- 
tion limit and the spectra suffer from strong contamination 
due to the diffuse emission. It is not clear whether it is the 
X-ray counterpart of the supernova remnant or not. 



5 DISCUSSION 

Long-term flux variability of individual sources could result 
from various kinds of physics and structure of the sources, 
revealing information about their physical nature. However, 
here we are concerned more about the source variability 
of the whole galaxy. Our study shows that there is an 
prominent fraction of transients and variables in M82. We 
are therefore interested in some questions: What kind of 
source population and mechanism could explain the long- 
term variability of M82? Is the source variability of star- 
forming galaxies higher than that of other types of galaxies? 
We believe that the overall variability of a galaxy reflects 
the nature of its underlying source population. Consider- 
ing that previous studies of X-ray sources in external galax- 
ies are hitherto mainly based on X-ray luminosity functions 
(XLFs), can we further constrain the source populations 
with the information of temporal properties? In this section, 
we will first discuss the source population and the possible 
origins of the long-term variability of M82. And then we 
will further compare the long-term variability of different 
morphological types of galaxies. 

Luminous X-ray sources in galaxies are believed to be 
dominated by XRBs. Based on previous studies of XLFs 
and other characteristics of our Galaxy and different types 
of external galaxies (see a review article by Fabbiano 2006), 
LMXBs are thought to dominate bright X-ray sources in 
early type galaxies. HMXBs dominate the source popula- 
tions of galaxies with high star formation rate (SFR) to to- 
tal stellar mass ratio. For spiral galaxies, the populations 
are more complex with a combination of both LMXBs and 
HMXBs, which differs from bulge, disk to arms. Providing 
the strong star-forming activities in M82, HMXBs associated 
with young stellar populations are expected to dominate the 
bright X-ray sources we detected. With XLFs of galaxies, it 
is possible to discriminate the contribution of HMXBs and 
LMXBs. The XLFs of HMXBs show a flatter slope of power 
law distribution without high luminosity cut-off. In contrast, 
high luminosity cut-off is seen in XLFs of LMXBs (Grimm 
et al. 2002). Grimm et al. (2003) further studied XLFs of 
star-forming galaxies. After excluding the estimated contri- 
bution from LMXBs, they suggested that there is a universal 
XLF of HMXBs characterized by a power law distribution 
with a cumulative slope of —0.61 ± 0.12. Here we present 
the XLFs of M82 in all the observations we used (Figure 4). 
Except for observation 6 and 7, the XLFs cover the sources 
within nearly all the -D25 isophote of M82. The low luminos- 
ity part of the XLFs of M82 is highly incomplete due to its 
strong diffuse emission. We set a completeness threshold at 
10''* ergs s~^and then flt the differential XLFs with power 
law using maximum likelihood. The cumulative slopes of all 
observations (except those with incomplete coverage of the 
D25 isophote) of M82 are shown in Table 3. The XLF of 
M82 shows a flat power law distribution without high lumi- 
nosity cut-off, suggesting a dominating HMXB population 
in the system. Furthermore, the good match of the cumula- 
tive slope and that of the universal XLF of HMXBs suggests 
that the contribution from LMXBs is fairly low. Therefore, 
we consider most of the sources in M82 to be HMXBs, which 
can constrain the origins of the long-term variability we de- 
tected. Apart from XLFs, we also consider a classification 
scheme developed by Prestwich et al. (2003) using X-ray col- 
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Figure 4. Cumulative X-ray luminosity functions of M82 in all 
observations. The observations cover nearly all the D25 isophotc 
of M82 except observation 6 and 7. 



Table 3. Cumulative slopes of the XLFs of M82 



No. 


Obs. ID 


Slope 


1 


361 


-0.72 


2 


1302 


-0.65 


3 


1411-1 


-0.64 


4 


1411-2 


-0.65 


5 


2933 


-0.65 


8 


8189 


-0.63 


9 


8505 


-0.63 



NOTE. — Observation 6 and 7 are not included because of the 
incomplete coverage of the D25 isophote of M82 . 

ors of sources. However, the high absorption of M82 makes 
it more difficult to probe the source population using X-ray 
colors. 

In addition, luminosity level reached by sources could 
impose constraint on source population since a higher ob- 
served luminosity of a source may imply a higher mass of its 
compact object if the Eddington limit is not violated. There- 
fore, luminosities above 10"^^ ergs s~^reduce the possibility 
for an XRB to host a neutron star. In fact, about half of 
the X-ray sources we detected in M82 exceed the Eddington 
limit of typical neutron stars (a few times 10^* ergs s~^), sug- 
gesting that there may be half or more BH XRB systems in 
M82. This statement can be further supported by the studies 
of source colors. Based on X-ray colors, Colbert et al. (2004) 
suggested that the dominant sources in nearby star-forming 
galaxies are generally inconsistent with accretion-powered 
NS HMXBs and more consistent with softer BH HMXBs. 
Besides, Swartz et al. (2004) showed that the properties of 
ULXs are similar to those of fainter sources, implying that 
ULXs are the high luminosity end of HMXBs, supernovae 
and LMXBs. It is therefore improper to treat ULXs as a 
distinct class of X-ray sources. In summary, we believe that 
the bright X-ray sources in M82 are dominated by HMXBs, 
and more than half of them should be black hole systems. 

Based on the above discussion about source population, 
we could further constrain the possible origins of the de- 
tected source variability into temporal phenomena related to 
HMXBs. Among different types of XRBs, various temporal 



properties have been observed and to some extent explained 
by theories. For Roche lobe overflow systems, BH LMXBs in 
the Galaxy are more variable than neutron star systems (see 
the X-ray variability functions of Galactic source population 
in Figure 5 of Feng and Kaaret 2006) . XRBs with a neutron 
star or a high mass stellar companion are thought to be rel- 
atively stable. Irradiation of the accretion disk resulted from 
the emission of a young massive donor star or the surface 
of a neutron star can maintain a temperature higher than 
the hydrogen ionization temperature, Th ^ 6500 K, thereby 
suppressing the instabilities in its accretion disk (e.g. van 
Paradijs 1996; Irwin 2006). On the other hand, HMXBs 
could also be transients. In our Galaxy, more than half of the 
HMXBs are systems consist of a Be star and a neutron star. 
The rapid rotation of the Be star created an extended cir- 
cumstellar disk. If the orbit is sufficiently eccentric, accretion 
could be triggered all of a sudden when the compact object 
pass through the disk. Therefore, almost all the Be HMXBs 
are transients. For HMXB dominant galaxies like M82, many 
transients should be Be HMXBs. Note that the compact ob- 
jects of Be HMXBs in our Galaxy are mainly neutron star, 
but more black hole systems are expected in M82. This does 
not oppose our inference. Considering the formation of Be 
HMXBs, the chance to form neutron star systems is fairly 
higher than that of black hole systems. Therefore, the pres- 
ence of some luminous black hole systems may imply a large 
amount of neutron star systems with lower luminosities. 

Based on strong source variability, we infer that many of 
the sources in star-forming galaxies are Be transients. This 
inference could be checked from other aspects. In a HMXB, 
mass is accreted through either stellar wind or Roche lobe 
overflow. For both the stellar wind systems and the Roche 
lobe overflow systems, the mass donor could be an OB- 
supergiant or a Be star. From a theoretical study (Postnov 
2003) of the universal XLF of HMXBs deduced by Grimm 
et al. (2003), the —0.6 cumulative slope of XLF could be 
explained by the mass-luminosity and mass-radius relations 
of wind-fed HMXBs. On the other hand, the slope of Roche 
lobe overflow systems is not consistent with the observed 
universal slope. Postnov (2003) considered that for HMXBs 
with Roche lobe overflow, the high efficiency supercritical ac- 
cretion produces luminosities on the order of the Eddington 
limit, which makes their XLF close to the mass distribution 
of the compact objects. The different XLF slope with re- 
spect to the —0.6 universal value suggests that Roche lobe 
overflow population may not make a significant contribu- 
tion. After limiting the major contribution to the univer- 
sal XLF in wind-fed HMXBs, we could further reject the 
wind-fed systems with a donor of supergiant, because they 
possess relatively lower luminosities (10"^^ — 10^^ ergs s~^) 
than the detected sources in our study. For these reasons, 
we exclude the significant contribution of Roche lobe over- 
flow systems and the fainter wind-fed supergiant systems. 
The major population of the X-ray sources in M82 and the 
observed universal HMXB XLF is most likely to be wind- 
fed Be HMXB. This is consistent with our previous inference 
based on source variability. 

Next, we discuss long-term source variability of galaxies 
with different morphological types. Not many studies about 
this have been done up to now. With our results of the proto- 
type star-forming galaxy M82 and previous studies of several 
galaxies, a rough picture could be constructed. 
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Figure 5. X-ray variability function of bright sources in M82 
combined with Feng and Kaaret's results (Figure 5 in Feng and 
Kaaret 2006) of the old elliptical galaxy NGC 1399 and the young, 
star-forming Antennae galaxies (NGC 4038/4039). 

The first way to compare source variability of differ- 
ent systems is based on X-ray variability function (XVF). 
Feng and Kaaret (2006) studied the temporal properties of 
bright X-ray sources in two extremely different host environ- 
ments: the old elhptical galaxy NGC 1399 and the young, 
star-forming Antennae galaxies (NGC 4038/4039). Maxi- 
mum ratios in luminosity, R, of sources were obtained using 
similar analysis process of our work. However, it should be 
noted that the luminosity ratios here are all lower limits. 
For both of their two target galaxies, they used 7 Chandra 
ACIS observations within 5 years. Here we present our X- 
ray variability function of M82 combined with the results 
of Feng and Kaaret's work (Figure 5). A luminosity cut-off 
of 3 X 10''* ergss"^ is set since the detection limits in most 
of the observations of these three galaxies are below this 
level. The early type galaxy, NGC 1399, hosts many bright 
X-ray sources, and their variability is clearly lower than the 
other two starburst systems. For early type galaxies, bright 
X-ray sources are dominated by LMXBs. Although many 
LMXBs in our Galaxy show transient phenomena, the bright 
sources in NGC 1399 are relatively stable with respect to 
those young HMXBs in star-forming galaxies. However, we 
expect that the variability of LMXBs could be more promi- 
nent on longer timescale. 

Compared with the Antennae galaxies, M82 hosts more 
X-ray sources beyond the threshold and more high R 
sources. However, it is premature to claim that the source 
variability of M82 is intrinsically higher than that of the 
Antennae. Because of the lower detection limit of M82, 
we can detect variations with higher R. It is possible that 
the detailed source populations and source variability are 
slightly different in star-forming galaxies with different con- 
ditions. To have a better understanding of the relation be- 
tween X-ray variability in star-forming environment and 
star- formation rate (SFR), star-formation history and/or 
other factors, more studies are required. 

We suggest another way to compare source variability 
of galaxies. As we mention in § 4.2, the fraction of variables 
with J? > 3 could be an indicator of variability. X-ray vari- 
ability functions plotted with the lower limit of R give a de- 
tailed picture of variability distribution. However, with cur- 
rent instruments, detection limits of external galaxies out- 
side the local group are not enough to cover the sources in 
their low states. Therefore, the XVFs, especially the high R 
part of them strongly depend on the detection limits. XVFs 



Table 4. Inferred properties of host galaxies 



Galaxy 


SFR" (A% yr-l ) 


Mass' (Mq) 


1VI82 


9.3 


6.3 X lO^O 


Antennae 


8.3 


5 X 10^° 


NGC 1399 




1.6 X 10^^ 


NGC 6946 


4 


2 X lO" 


NGC 4485/4490 


4.7 


1.8 X 10^" 



" M82, Antennae (Colbert et al. 2004), NGC 6946 (Sauty et 
al. 1998) and NGC 4490 (Clemens et al. 1999). 
*' Stellar mass of M82, Antennae and NGC 1399 were cal- 
culated using Lj; (Colbert et al. 2004). Dynamical mass is 
used in NGC 6946 (Crosthwaite and Turner 2007) and NGC 
4485/4490 (Clemens et al. 1999). 

should only be used to compare systems with similar detec- 
tion limits. To compare the few studies done with different 
detection limits, the fractions of variables with i? > 3 pro- 
vide a rough but proper description. 

Two systems with moderate star-forming activities, 
NGC 6946 and NGC 4485/4490, were well studied with 
an emphasis on long-term variability by Fridriksson et al. 
(2008). The detailed data analysis and the detection limits 
(a few times lO^'' ergss"^ for NGC 6946 and ~ 10^^ ergss"^ 
for NGC 4485/4490) are not exactly the same as ours but 
comparable. We did not compare their results with ours us- 
ing XVF since there are relatively few bright sources beyond 
the threshold we set in the previous comparison. We recal- 
culated the fractions of variables with i? > 3 of these two 
systems using the criteria we used in this paper. We ob- 
tained 23% for NGC 6946 and 21% for NGC 4485/4490. 
The long-term variability of the X-ray sources in M82 (47% 
sources are variables with flux ratio > 3) is extremely high 
compared to these two galaxies with star-forming activities. 
Since the SFRs in these two galaxies are lower than that of 
M82 (see the SFRs of galaxies in Table 4), it is reasonable to 
claim that SFR may be a correlated factor of long-term X- 
ray variability. This is also consistent with the concern that 
there is a causality between source populations and tempo- 
ral properties since SFR is an important indicator of source 
populations. 



6 CONCLUSIONS 

We have reported a study of the long-term variability of the 
X-ray sources in the starburst galaxy M82. By combining 
9 archival Chandra data sets observed between 1999 and 
2007, we detected 58 X-ray point sources within the -D25 
isophote of M82 down to a luminosity of ~ 10"^^ ergs s~^. 
Seven detections can be expected from background sources. 
Three sources are coincident with the positions of historical 
supernovae or radio-selected SNRs. 

Twenty-six sources in M82 exhibit long-term (i.e., days 
to years) flux variability and 3 exhibit long-term spectral 
variability. Furthermore, we classified 26 sources as variables 
and 10 as persistent sources. Among the total 26 variables, 
17 varied by a flux ratio of > 3 and 6 are transient candi- 
dates. To reduce the selection effects, we used the fraction 
of variables with a flux ratio of > 3 to be an indicator of 
long-term variability. The long-term variability of the X-ray 
sources in M82 is higher than Early-type galaxies and some 
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galcixies with star-forming activities (NGC 6946 and NGC 
4485/4490). 

For the variables, transient candidates and spectral vari- 
ables we classified, various kinds of variability properties are 
shown. Two of the transients showed recurrent outbursts. 
One of the variables was classified as SSS. M82 X-f showed 
a transition from a low/hard state to a high/soft state in 
its color-luminosity diagram, and another spectral variable 
showed an inverse correlation. 

Based on the strong variability, the X-ray luminosity 
functions, and the number of ULXs, we suggest that the X- 
ray source population of M82 with luminosities greater than 
~ 10^^ ergs s"^ is dominated by HMXBs. More black hole 
XRBs should be expected than neutron star XRBs. Some of 
the transients should be Be HMXBs. 
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